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Summary. — 

Chandrasekhar separated the Dirac equation for spinning and massive particles in 
Kerr geometry into radial and angular parts. In the present review, we present 
solutions of the complete wave equation and discuss how the Dirac wave scatters off 
Kerr black holes. The eigenfunctions, eigenvalues and reflection and transmission co- 
efficients are computed for different Kerr parameters. We compare the solutions with 
several parameters to show how a spinning black hole distinguishes mass and energy 
of incoming waves. Very close to the horizon, the solutions become independent of 
the particle parameters indicating an universality of the behaviour. 

PACS 04.20.-q - Classical general relativity. 
PACS 04.70.-S - Physics of black holes. 
PACS 04.70.Dy - Quantum aspects of black holes. 
PACS 95.30. Sf - Relativity and gravitation. 



II Nuovo Cimento (in press) 
1. Introduction 



One of the most important solutions of Einstein's equation is that of the spacetime 
around and inside an isolated black hole. The spacetime at a large distance is flat 
and Minkowskian where usual quantum mechanics is applicable, while the spacetime 
closer to the singularity is so curved that no satisfactory quantum field theory could be 
developed as yet. An intermediate situation arises when a weak perturbation (due to, say, 
gravitational, electromagnetic or Dirac waves) originating from infinity interacts with a 
black hole. The resulting wave is partially transmitted into the black hole through the 
horizon and partially scatters off to infinity. In the linearized ('test field') approximation 
this problem has been attacked in the past by several authors [1, 2, 3, 4]. The master 
equations of Teukolsky [2] which govern these linear perturbations for integral spin (e.g., 
gravitational and electromagnetic) fields were solved numerically by Press & Teukolsky 
[5] and Teukolsky & Press [6] . While the equations governing the massive Dirac particles 
were separated in 1976 [3], the angular eigenfunction and eigenvalue (which happens 
to be the separation constant) have been obtained in 1984 [7] and radial solutions have 
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been obtained only recently [8, 9, 10, 11]. Particularly interesting is the fact that whereas 
gravitational and electromagnetic radiations were found to be amplified in some range 
of incoming frequencies, Chandrasckhar [4] predicted that no such amplifications should 
take place for Dirac waves because of the very nature of the potential experienced by the 
incoming fields. Although this later conclusion was drawn using an asymptotic equation, 
wc show that this is indeed the case even when complete sohitions arc considered for the 
Dirac wave perturbations. Chandrasckhar also speculated that one needs to look into 
the problem for negative eigenvalues {X) where one might come across super-radiance for 
Dirac waves. 

In the present review, we discuss this important problem and its solutions. We show 
the nature of the radial wave functions as a function of the Kerr parameter, rest mass 
and frequency of incoming particle. We also verify that super-radiance is indeed absent 
for the Dirac field. Unlike earlier works [5, 6] where numerical (shooting) methods were 
used to solve the master equations governing gravitational and electromagnetic waves, 
we use a classical approach whereby we approximate the potential felt by the particle by 
a collection of small steps. 

Below, we present the separated Dirac equations from Chandrasckhar [4] using the 



same choice of units: we choose H = 1 = G = c, so that the unit of mass becomes 



The equations governing the radial wave- functions Rj.i corresponding to spin ±5 
respectively are given by: 




the unit of time becomes 




and the unit of length becomes 




AiDoi?_i 



1 = {^ + impr)AiR_^_i, 



(la) 



AivlAiR_^_i = (:X - impr)R_i, 
where, the operators Vn and Vj^ are given by. 



(16) 



1^n = dr + -^+ 2n 



(r-M) 
A 



(2a) 



{r-M) 
A 



(26) 



and 



A = + a^ - 2Mr, 



(3a) 



K = {r"^ + a'^)a + am. 



(36) 



Here, a is the Kerr parameter, n is an integer or half integer, a is the frequency of incident 
wave, M is the mass of the black hole, nip is the rest mass of the Dirac particle, % is the 
eigenvalue of the Dirac equation and m is the azimuthal quantum number. 
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The equations governing the angular wave- functions S^i corresponding to spin ±i 



respectively are given by: 



j0.iS_^_i = — {X— amp cos 9)S_ 



(4a) 



jC\S_i = +{'X + ampCos9)S,i 
where, the operators and £^ are given by, 



(46) 



Cn = de + Q + ncot 6, 



(5a) 



and 



Cl^ = d$ — Q + ncot( 



Q = aasinO + m cosec 6. 



(56) 



(6) 



Note that both the radial and the angular sets of equations i.e., eqs. l(a-b) and eqs. 
4(a-b) are coupled equations. Combining eqs. 4(a-b), one obtains the angular eigenvalue 
equations for the spin-^ particles as [7] 



mil fi tiiiiO ^-j- 



C\ + {X^ - a^ml cos^ ■ 



23 ^ + amp cos 6 2 



S_i = 0. 



(7) 



There are exact solutions of this equation for the eigenvalues X and the eigenfunctions 
S_i when p = ^ = 1 in terms of the orbital quantum number I and azimuthal quantum 
number m. These solutions are [7]: 



= (Z + + aa{p + 2m) + a^a^ 



1 - 



y 



2{l + 1) + aax 



and 



where. 



aay 



2{l + 1) + aax 5 



(8) 



(9) 



p = F{l,l); x^F{l + 1,1 + 1); y = F{l,l+l) 



and 



F{h,l2) = [(2/2 + l)(2Zi + !)]*< l2'^mO\km > 
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(10) 

with < ....|.. > are tlie usual Clebsh-Gordon coefficients. For other values of p one has 

to use perturbation theories. Sohitions upto sixth order using perturbation parameter 
aa are given in Chakrabarti [7] and are not described here. The eigenfunctions X are 
required to solve the radial equations which we do now. 

The radial equations 1 (a-b) are in coupled form. One can decouple them and express 
the equation either in terms of spin up or spin down wave functions R±^ but the expres- 
sion loses its transparency. It is thus advisable to use the approach of Chandrasekhar [4] 
by changing the basis and independent variable r to, 



2Mr+ + am/a^ 

r* = r H log 

r+ — r_ 



r 



2Mr- + am/a 
r+ — r_ 



log 



1 (r>r+). (11) 



where, 



_d_ 

dr* 



2 I 2 



+ am/a, 



(12) 



to transform the set of coupled equations l(a-b) into two independent one dimensional 
wave equations given by: 



A* 



— -^a]p^^= —(X - impr)P_i; [—+ia]P_i= —{X + impr)P+i 



(13) 

Here, Vq = ^{-£- + ia) and Vq = ^{-^ ia) were used and wave functions were 

redefined as R_i — P_i and AiR.i = P.i. 

2 2 ''2 "'"2 

We are now defining a new variable, 

e = tan~^{mpr/X) (14) 

which yields. 



cosf = 



V(X2 + m2r2)^ 



and sin^ = 



Vi^"" + mlr^) 



and 



{X ± impv) = exp{±ie)^{:>? + mlr^), 



so the coupled equations take the form 

Ai 



_d_ 



^ 9 



{X^ +mlr^)^/^P_iexp 



-i tan 



mpT 



(15a) 



and 



A* 



-i-+ia)P_i = -^{X^ + my)^/^P:iexp 



itan 



(156) 



DIRAC EQUATION IN KERR GEOMETRY AND ITS SOLUTION 



5 



Then defining, 



and 

we obtain, 



dr* 



P- 1 = V'- i exp 



. / A "Xrrip 1 
^c^ 1 + 



1 . , _i / mpr 
tan — — 
2 \ X 



1 ; iii'pi 



-—I tan 
2 V X 



(16a) 



(166) 



w2 2a ^^2+^2^2 



^+1 =-^(^'+my)i/2V-i (17a) 



and 



+ ia I 1 + 



dr* I a;2 2a 7? + mlr^j = ^^2 v p / r+. 



Further choosing f* = r* + ^tan i(2^) so that (ir* = (1 + ;^^^^^t^:^^)c^''*) the 
above equations become. 



(176) 



and 



where, 



(18a) 



(186) 



W = 



uj'^iX^ + mlr'^) + XmpA/2a ' 
Now letting Z± = V'+i ± V'-i we can combine the differential equations to give 



(19) 



— -W\Z+^ iaZ_, 
dr* 



(20a) 



and 



+ W] Z_ = iaZ+. (206) 

From these equations, we readily obtain a pair of independent one-dimensional wave 
equations. 
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where, V± = ± 



Ai {^^ + my )^/^ 



a'\Z± = V±Z±. (21) 



[w^{X^ + my) +XmpA/2a]- 



:[Ai{X'' + myf^ ± ((r - M){7? + my) + Sm^rA)] 



2r^2^ A/9 13 P^(^' + "^^') + + - M)l<y\. (22) 

One important point to note: the transformation of spatial co-ordinate r to r, (and f*) 
is taken not only for mathematical simplicity but also for a physical significance. When 
r is chosen as the radial co-ordinate, the decoupled equations for independent waves 
show diverging behaviour. However, by transforming those in terms of r* (and f*) we 
obtain well behaved functions. The horizon is shifted from r = r_|_ to f* = — oo unless 
a < as = —am/2Mr+ (eq. 11). In this connection, it is customary to define CTc where 
= (eq. 13). Thus, CTc = —m/a. If cr < a^, the super-radiation is expected [4]. 

2. — Solution Procedure 

The choice of parameters is generally made in such a way that there is a significant 
interaction between the particle and the black hole, i.e., when the Compton wavelength 
of the incoming wave is of the same order as the outer horizon of the Kerr black hole. 
Similarly, the frequency of the incoming particle (or wave) should be of the same order 
as inverse of light crossing time of the radius of the black hole. These yield [8], 



m„ 



a [M + ^/lM^ - a'^)]-\ (23) 



One can easily check from equation (22) that for r — > oo (i.e., f* — > oo) V± — > rrip. So 
the total energy of the physical particle should greater than square of its rest mass. So 
if we expand the total parameter space in terms of frequency of the particle (or wave), 
a and rest mass of the particle, rup, it is clear that 50% of total parameter space where 
(T < rrip is unphysical (In this case, the energy is such that a particle released from a 
finite distance cannot go back to infinity after scattering.), and one need not study this 
region. Out of the total physical parameter space there are two cases of interest: (1) the 
waves do not 'hit' the potential barrier and (2) the waves do hit the potential barrier. 
To sole these potential problem first, we replace the potential barrier by a large number 
of steps as in the step-barrier problem in quantum mechanics. Fig. 1 shows one such 
example of the potential barrier [10] V+ (Eq. 22) which is drawn for a = 0.5, mp = 0.8 
and a = 0.8. In reality we use tens of thousands of steps with suitable variable widths 
so that the steps become indistinguishable from the actual function. The solution of Eq. 
22 at nth step can be written as [14], 

Z+^n = Anexp[iknf*,n] + Bnexp[-iknh,n] (24) 
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Fig. 1 : Behaviour of V+ (smooth solid curve) for a = 0.5, rUp — 0.8, a = 
0.8. This is approximated as a coUection of steps. In reality tens of thousand 
steps were used with varying step size which mimic the potential with arbitrary 
accuracy. 



when energy of the wave is greater than the height of the potential barrier. The standard 
junction condition is given as [14], 

Z+,n = Z+ n+1 and = —rr-\n+l- (25) 



The reflection and transmission co-efficicnts at nth junction are given by: 

An+l{kn+i - kn) + Bn+i{kn+i + kn) , 
An+l{kn+l + kn) + iJ„+l(A:„+l - kn) 

At each of the n steps these conditions were used to connect solutions at successive steps. 
Here, k is the wave number (fc = y'<J^-V±) of the wave and kn is its value at nth step. 
We use the 'no-reflection' inner boundary condition: i? ^ at r» — > —oo. 
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Fig. 2a: Reflection (R) and transmission (T) coeflicients of waves with varying 
mass as functions of r«. rup — 0.78 (solid), rup — 0.79 (dotted) and mp = 0.80 
(long-dashed) are used. Other parameters are a — 0.5 and a — 0.8. Inset shows 
R in logarithmic scale which falls off exponentially just outside the horizon. 



For the cases where waves hit on the potential barrier, inside the barrier (where 
(T^ < V+) we use the wave function of the form 

Z+_„ = Anexp[-anf^^n] + Bnexp[anf*,„] (27) 



where, q;„ — \/V± — cr^, as in usual quantum mechanics. 

3. Examples of Solutions 

Fig. 2a shows three solutions [amplitudes of Re(2'_(.)] for parameters: a = 0.5, a = 0.8 
and mp — 0.78, 0.79, and 0.80 respectively in solid, dotted and long-dashed curves. The 
energy is always higher compared to the height of the potential barrier (Fig. 1) and 
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Fig. 2b: Amplitude of Re(Z+) of waves with varying mass as functions of f«. 
rup — 0.78 (solid), nip — 0.79 (dotted) and mp = 0.80 (long-dashed) are used. 
Other parameters are a — 0.5 and a — 0.8. 



therefore the particles do not 'hit' the barrier, k goes up and therefore the wavelength 
goes down monotonically as the wave approaches a black hole. It is to be noted that 
though ours is apparently a 'crude' method, it has flexibility and is capable of presenting 
insight into the problem, surpassing any other method such as ODE solver packages. 
This is because one can choose (a) variable steps depending on steepness of the potential 
to ensure uniform accuracy, and at the same time (b) virtually infinite number of steps 
to follow the potential as closely as possible. For instance, in the inset, we show R in 
logarithmic scale very close to the horizon. All the three curves merge, indicating that 
the solutions are independent of the mass of the particle and a closer inspection shows 
that here, the slope of the curve depends only on a. The exponential dependence of i?„ 
close to the horizon becomes obvious. Asymptotically, V± = nip (eq. 22), thus, as nip 
goes down, the wavelength goes down. In Fig. 2b, we present the instantaneous values 
of the reflection R and transmission T coefficients (i.e., i?„ and r„ of Eq. 26) for the 
same three cases. As the particle mass is decreased, k goes up and corresponding R goes 
down consistent with the limit that as fc — *■ oo, there would be no reflection at all as in 
a quantum mechanical problem. 



10 



SANDIP K. CHAKRABARTI AND BANIBRATA MUKHOPADHYAY 




Fig. 3a: Reflection (i?) coefficient of waves witli varying mass as functions 
of f*. nip = 0.16 (solid), nip = 0.164 (dotted) and nip — 0.168 (long- 
dashed) are used. Other parameters are a = 0.95 and a = 0.168. 



Figs. 3(a-b) compare a few solutions where the incoming particles 'hit' the potential 
barrier. We choose, a = 0.95, a = 0.168 and mass of the particle nip = 0.16, 0.164, 0.168 
respectively in solid, dotted and long-dashed curves. Inside the barrier, the wave decays 
before coming back to a sinusoidal behaviour, before entering into a black hole. In Fig. 
3b, we plotted the potential (shifted by 2.05 along vertical axis for clarity). Here too, 
the reflection coefficient goes down as k goes up consistent with the classical result that 
as the barrier height goes up more and more, reflection is taking place strongly. Note 
however, that the reflection is close to a hundred percent. Tunneling causes only a few 
percent to be lost into the black hole. 

Figs. 4(a-b) show the nature of the complete wave function when both the radial and 
the angular solutions [7] are included. Fig. 4a shows contours of constant amplitude 
of the wave (i?_i/2»S'-i/2) in the meridional plane - X is along radial direction in the 
equatorial plane and Y is along the vertical direction. The parameters are a = 0.5, 
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Fig. 3b: Amplitude of Re{Z+) of waves with varying mass as functions of 
f*. irip = 0.16 (solid), mp — 0.164 (dotted) and nip — 0.168 (long-dashed) 
are used. Nature of potential with mp = 0.168 is drawn shifting vertically 
by 2.05 unit for clarity. Other parameters are a — 0.95 and cr = 0.168. 



mp — 0.8 and a = 0.8. Some levels are marked. Two successive contours have amplitude 
difference of 0.1. In Fig. 4b a three-dimensional nature of the complete solution is given. 
Both of these figures clearly show how the wavelength varies with distance. Amplitude 
of the spherical wave coming from a large distance also gets weaker along the vertical 
axis and the wave is forced to fall generally along the equatorial plane, possibly due to 
the dragging of the inertial frame. 

4. Conclusion 

We review here the scattering of massive, spin-half particles from a spinning black 
hole with particular emphasis to the nature of the radial wave functions and reflection 
and transmission coefficients. Here we presented a well known quantum mechanical step- 
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Fig. 4a: Contours of constant amplitude are plotted in the meridional plane 
around a black hole. Radial direction on equatorial plane is along X axis and 
the vertical direction and along Y. Both radial and theta solutions have been 
combined. Parameters are a = 0.5, nip = 0.8 and a = 0.8. 



potential approach [10] but one can verify by any numerical technique that the solution 
would remain the same. A modified WKB approximation [8, 9, 11] also yields similar 
result in Kerr geometry [15]. The approach presented here (i.e., step potential approach) 
is very transparent since a complex problem of barrier penetration in a spacetime around 
a spinning black hole could be tackled very easily. We report a few significant observa- 
tions of these papers that the wave function and R, and T behave similarly close to the 
horizon independent of the initial parameter, such as the particle mass rup. Particles of 
different mass scatter off to a large distance completely differently, thus giving an im- 
pression that a black hole could be treated as a mass spectrograph! When the energy of 
the particle becomes higher compared to the rest mass, the reflection coefficient dimin- 
ishes as it should it. Similar to a barrier penetration problem, the reflection coefficient 
becomes close to a hundred percent when the wave hits the potential barrier. Another 
significant observation is that the reflection and transmission coefficients are functions 
of the radial coordinates. This is understood easily because of the very nature of the 
potential barrier which is strongly space dependent which we approximate as a collection 
of steps. Combining with the solution of theta-equation, we find that the wave-amplitude 
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Fig. 4b: Three dimensional view of Ji_ 1/2 5_ 1/2 are plotted in the meridional 
plane around a black hole. Both radial and theta solutions have been combined. 
Parameters are a = 0.5, nip = 0.8 and a = 0.8. 



vanishes close to the vertical axis, possibly due to the frame-dragging effects. 
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